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1 Introduction 

In a number of recently published papers [H [9l [TOl [19] we have developed a distinctive geometric 
approach to the study of regular Lagrangian dynamical systems, and especially to the problem 
of formulating reduced equations for systems which are invariant under the action of a symmetry 
group. The main distinctive features of our approach are, firstly, the formulation of the Euler- 
Lagrange equations in a way which is well adapted to the idea that their function is to determine 
a vector field on the velocity phase space which is of second-order type (so that the differential 
equations which determine its integral curves are of second order in the configuration space 
coordinates), and yet is completely coordinate independent; and secondly, the consistent use of 
anholonomic frames and their associated quasi- velocities. In this paper we shall extend these 
ideas to cover Lagrangian systems subject to nonholonomic linear constraints, for which both 
the Lagrangian function and the constraint distribution are invariant. Such constraints arise 
naturally in the context of systems with rigid bodies rolling without slipping over a surface or 
possessing a contact point with the surface in the form of a knife edge. A classical reference for 
the dynamics of mechanical systems with nonholonomic constraints is the book by Neimark and 
Fufaev [21] . The recent books [H [71 [12] contain many references to different modern geometric 
approaches to the theory. We will work with autonomous systems; for formulations of the 
nonholonomic dynamics in a time-dependent set-up see e.g. \16\ 122], 

The formulation of the Euler-Lagrange equations mentioned above goes as follows. We consider 
a Lagrangian system over a differentiable manifold Q (configuration space). The Lagrangian L 
is a function on the tangent bundle r : TQ — t- Q (velocity phase space); it is regular if its Hessian 
with respect to the fibre coordinates is nonsingular. The following proposition holds jlOj . 



Proposition 1. Let L be a regular Lagrangian on TQ. There is a unique second-order differ- 
ential equation field T such that 

r(Z^(L)) - ZC(L) = 
for all vector fields Z on Q. Moreover, V may he determined from the equations 

r(X7(L)) - Xf{L) = 0, i = 1, 2, . . . n = dimQ, 
for any frame {Xi} on Q (which may he a coordinate frame or may he anholonomic) . 

Here and Z'^ are respectively the vertical lift and the complete or tangent lift of Z to TQ 
(we refer to [iT] for the most common notions on tangent bundle geometry). The formulation 
in terms of the frame {Xi} leads directly to Hamel's equations 

where Xi = X^d/dq^; the are the quasi- velocities associated with the frame; and the coeffi- 
cients R^j are defined by [Xi,Xj] = RfjX^ and are collectively called the object of anholonomity 
of the frame. We refer the reader to Section |4] for more details. 

The equations determining the dynamics of a regular system subject to nonholonomic linear 
constraints admit a rather similar formulation. The constraints may be specified in either of two 
equivalent ways: as a distribution D on Q (the constraint distribution), or as a submanifold C of 
TQ (the constraint submanifold). The two are related as follows: C = {(g, u) G TQ : u £ T>q C 
TqQ}. We assume that the dimension of each Vg, and equivalently the fibre dimension of Cg, is 
constant and equal to m. A vector field T on C is said to be of second-order type if it satisfies 
T*{q,u)'^ = u for all {q, u) £ C. A Lagrangian function L is said to be regular with respect to T> 
if for any local basis {^a} of V, 1 < a < m, the symmetric m x m matrix whose entries are 
X^{XJ{L)) (functions on C) is nonsingular. In [TO] we proved the following proposition. 

Proposition 2. Let L he a Lagrangian on TQ which is regular with respect to T). Then there is 
a unique vector field T onC which is of second- order type, is tangent to C, and is such that on C 

T{Z^{L))- Z''{L) = {) 

for all Z £T>. Moreover, T may he determined from the equations 

r(X^(L))-XS(L) = 0, a = l,2,...m, 

on C, where {X^} is any local hasis for T>. 

This is our version of the Lagrange-d'Alembert principle (see [HE] for other versions); the vector 
field r is the dynamical field of the constrained system. 

The formal similarity between the standard Euler-Lagrange equations and the Lagrange-d'Alembert 
equations in these formulations is self-evident. We shall exploit this similarity in deriving the 
reduced equations for an invariant constrained system: as we shall show, to obtain those equa- 
tions it is enough to follow the reduction procedure for an invariant unconstrained system, while 
restricting attention to the constraint submanifold. 
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There are in fact two well-known ways of reducing the equations of an invariant unconstrained 
Lagrangian system. One method does not even take the Lagrangian structure of the system 
into account and simply involves factoring out by the action of the group, and leads to the 
so-called Lagrange-Poincare equations; this is described in e.g. [5l 1131 [T9]. The second does take 
advantage of momentum conservation; its first step is to restrict to a level set of momentum, 
and this is followed by a reduction with respect to the invariance group of the chosen value of 
the momentum. This is a generalized version of Routh's procedure; it is discussed in [8l[l8]. For 
a constrained system, however, these two methods are not equally applicable. This is because, 
although the constraint distribution is invariant under the symmetry group, it is not usually 
the case that any fundamental vector field of the action belongs to it. There is consequently 
no conservation of momentum, and no possiblility of Routh-type reduction. The greater part of 
this paper is therefore devoted to the adaptation of Lagrange-Poincare reduction to constrained 
systems. This we discuss in full generality: whereas many other papers ( [H [3l [6l [20] for example) 
restrict their attention to the case in which at each point q Q the constraint distribution Vq 
and the tangent space Vq to the orbit of the action together span TqQ, we make no such so-called 
'dimension assumption'; our only requirement is that Pg n Vq has constant dimension. 

Though Routh-type reduction is not possible in general, it can arise in particular cases, where 
there is a Lie subgroup H of the symmetry group G, necessarily normal, with Lie algebra f), such 
that for all ^ G {) the corresponding fundamental vector field ^ lies in V. Symmetries belonging 
to H are said to be horizontal. We devote a separate section to the discussion of this case. 

The paper is laid out as follows. In the following section we deal with the fundamental definitions 
and results concerning invariance of a constrained system under the free and proper action of a 
Lie group G on Q, leading to a version of the Atiyah sequence for such a system. In Section 3 
we give a resume in general terms of the Lagrange-Poincare reduction procedure, and show how 
it may be adapted to the case of an invariant constrained system. In Section 4 we derive explicit 
Hamel-type formulae, in terms of a (possibly) anholonomic frame, for the Lagrange-Poincare 
equations and the Lagrange-d'Alembert-Poincare equations successively. In Section 5 we discuss 
Routh-type reduction for systems with horizontal symmetries. 

2 Invariance of nonholonomic systems 

Assume from now on that a connected Lie g roup G cicts in a free and proper way on the left on 
the configuration manifold Q. Then vr : Q — )■ Q/G is a principal fibre bundle. The action ipg on 
Q induces an action Tipg on TQ. We will write A for the infinitesimal generator of the action 
on Q, associated to a Lie algebra element ^4 € g. Then A'~^ is an infinitesimal generator for the 
action on TQ. As in e.g. [3], we say that the nonholonomic system is invariant under G, or that 
it admits G as a symmetry group, if both the Lagrangian L and the constraint submanifold C 
of the system are invariant under the induced action of G on TQ. 

Proposition 3. The constraint submanifold C C TQ is invariant under Tip if and only if the 
constraint distribution T> on Q is invariant under ip. 

Proof. Since C = {{q, u) : u (z T>q}, C is invariant under Tip if and only if for every q Q, u & T>q 
and g & G, ipg^fU € V^^f^qy, that is to say, for every q G Q and g £ G, T'^gi^q) = ipg^Vq. □ 

Proposition 4. If L is regular with respect to T> the vector field T is invariant under the induced 
action of G on C. 
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Proof. For any A € and Z ^ D we have 

= i^(r(Z^(L))-Z°(L)) 

= [i^r](z-(L)) -r(ic(z-(L))) - [A^,Z%L) 
= [i^r](z-(L)) -r([i,z]-(L)) - lAzfiL). 

Now [A, Z] € V, due to the assumed hivariance of D. By the Lagrange-d'Alembert equation the 
last two terms above vanish. On the other hand, the bracket [A'-^, F] is vertical. This is certainly 
true for a second-order differential equation field on TQ, by a simple calculation in coordinates. 
Now r is a second-order differential equation field on C; but we can evidently extend it to a 
second-order differential equation field on a neighbourhood of C in TQ. Since both A'^ and T are 
tangent to C, so also is their bracket. So on C, [A'^,T] is independent of the choice of extension, 
and is vertical. It follows from the fact that the equation r](Z^(L)) = holds for all Z G P, 
and the assumption that L is regular with respect to T>, that [74'^,r] = 0. This may easily be 
seen by expressing everything in terms of the vertical lifts of a local basis for V. But [^4'-', T] = 
is the infinitesimal condition for T to be invariant. □ 

Since T is invariant, it reduces to a vector field T on C/G. Our main overall aim in this paper 
is to show how to determine this reduced vector field. We begin however by considering some 
aspects of the structure of nonholonomic systems which are invariant in the sense defined above. 

As we have already noted, tt : Q Q/G is a principal fibre bundle. Since T> is invariant it 
defines a distribution V on Q/G hy '^niq) = '^*(J^q)] this is well-defined because tt o ipg = n. Let 
us assume that T> has constant dimension. Then Q/G is equipped with a regular distribution 
P. Denote the corresponding submanifold (indeed vector subbundle) of T{Q/G) by C. 

Let Vq = ker7r^,g. Note that kervr^^^jx)^ = Vg n Vq. Let us denote it by Sg. Evidently S is an 
invariant distribution on Q, which is of constant dimension by the corresponding assumption 
for v. Since Sg C Vg for each g € Q, we may identify Sg with a vector subspace g'' of g, where 
q'^ = {A £ q \ Ag ^ Sg}. In terms of TQ, we can express g"? as follows. For w £ TgQ, w £ Vg 
if and only if is tangent to C; thus A £ q'^ and only if Ag is tangent to C. Since (see e.g. 

El El) 

we have 

qMi) = ad(g-i)g^. 

Consider = {{q,A) \A € g*^}. There is an action of G on g^ given by 

{q,A)^{i;g{q),ad{g~')A). 

On taking the quotient we obtain a vector bundle over Q/G, say , which is a vector subbundle 
of g = (Q X q)/G Q/G, the adjoint bundle associated with the principal G-bundle Q. 

Proposition 5. We have the following short exact sequence of vector bundles over Q/G: 

O^f^^C/G^C^O. 

This is a version for constrained systems of the so-called Atiyah sequence (see e.g. [9l I13|). 

^ g ^ TQ/G ^ T{Q/G) ^ 0. 

Each term in the sequence of the proposition is a subbundle of the corresponding term in the 
Atiyah sequence. In the next section we will use this observation when we divide the reduced 
equations for an invariant nonholonomic system into two sets. 
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3 Lagrange-Poincare reduction: generalities 



3.1 Standard Lagrange-Poincare reduction 

Before considering reduction of invariant nonholonomic systems we discuss Lagrange-Poincare 
reduction of the standard Euler-Lagrange equations. Recall the Euler-Lagrange equations as 
they appear in Proposition [H Assume that L is G-invariant: then so is P; it reduces to a vector 
field r on TQ/G, which we want equations for — so-called reduced equations. 

There is a sense in which the reduction of the Euler-Lagrange equations is immediate. The key 
step is to rewrite them in G-invariant form. It is enough to take Z to be invariant. Then 
and Z'^ are invariant under the induced action of G on TQ, and so define vector fields on TQ/G^ 
which we denote by Z^ and Z'^, though of course they are not vertical or complete lifts. The 
function Z'-^'{L) is invariant, and so defines a function on TQ/G^ which is just Z'^{1) (where I 
is the reduced function of L on TQ/G); likewise for Z^{L). Then the reduced equations are 
simply 

f - = 0, 

on TQ/G, for all invariant vector fields Z on Q] they are called the Lagrange-Poincare equations. 

However, we can be more explicit. The Euler-Lagrange equations can be divided into two sets, 
according to whether we take Z to be tangent to the fibres of vr or transverse to them. 

The Lagrange-Poincare equation for momentum. Consider first the Euler-Lagrange equa- 
tion T[Z^{L)) — Z'^(L) = for any vector field Z which is vertical with respect to n : Q —?■ Q/G. 
Such a vector field is determined by a g-valued function on Q, where Zg = C{Q)q- The mo- 
mentum p is a g*-valued function on TQ, which is G-equivariant under the usual action of G on 
TQ and the coadjoint action on q* . 

Take first A ^ g. We have 

(as real-valued functions on TQ; the angle brackets denote the pairing of q and q*). The 
conservation of momentum is just T{A,p) = (the Euler-Lagrange equation with Z = A): or, 
since A is constant and arbitrary, T{p) = 0. 

Now consider {C,p)- applying Leibniz' rule we have 

r(C,p) = (r(C),p) + (C,r(p)) = (c,p) 

(using the fact that r(/) = / for any function / on Q). We claim that the (almost tautological) 
equation T{(^,p) = {C,p) is the Euler-Lagrange equation for vertical Z. We compute Z^ and Z'^ 
in terms of C, as follows. Take a basis {Er} of g, and set (" = C^Er, where the coefficients C*" are 
functions on Q. Then Z = C^Er, and so 

= CE^, Z^ = CE^ + CE^- 

Then 

Z-(L) = CPr = (Cp), Z%L) = CPr = {C,P) 

as claimed. 

For reduction we need to take Z to be G-invariant: it will be so if and only if (" is G-equivariant 
(with now the adjoint action on g); this amounts to taking a section of the adjoint bundle 
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— >■ Q/G over Q. It is then clear that (CjP) wiU be invariant. It must then be the case (from 
the Euler-Lagrange equation and the invariance of F) that {C,p) is invariant. This can be shown 
directly too, in various ways. Here is a vector field version. The G-equivariance of p can be 
expressed as 

{B,A%p)) = -{[A,Blp), 

where A^B ^ q and i?] is their bracket in g. The assumed G-equivariance of Q is just 
A{C,) = [A,Q\ (again, bracket in q). Then 

A^tp) = {A%0,p) + {C,A%p)) 

= ^|(i(C)),p^-([AC],p) 
= ^|([AC]),p)-([AC],P> = o, 

using the obvious fact that X'^{f) = d/dt{X{f)). 

Since {C,p) and {C,p) are invariant they define functions on TQ/G, which we denote by {{C,p)) 
and {{C,p))- The corresponding reduced equation is 

r{{c,p)) = {{C,p)). 

We call it the Lagrange-Poincare equation for momentum. 

The horizontal Lagrange-Poincare equation. To obtain an invariant Euler-Lagrange equa- 
tion corresponding to the transverse directions we can use a principal connection on vr : Q — )> 
Q/G, or, equivalently, a splitting of the Atiyah sequence. Suppose we have such a connection: 
for any vector field Y on Q/G let be its horizontal lift to Q; then the transverse (let's call 
it horizontal) Euler-Lagrange equation is 

r((y«)^(L)) - (y«)^(L) = o. 

Incidentally, this expression is C°°((5/G)-linear in Y. Moreover, each term is G-invariant. 

We will express this equation in a different way. First we recall the construction of the Vilms 
connection from [9] (which is a special case of a more general construction in |23] ) . The complete 
lift of a type (1,1) tensor T on Q is given by |11] 

r^(x^) = T{xf, T°(x^) = T{xy. 

The original connection on vr : Q — )> Q/G can be represented by a type (1, 1) tensor field uj, so 
that uj{X) = if and only if X is horizontal, u}{V) = V for V vertical. Then the type (1, 1) 
tensor on TQ defining the Vilms connection is just cj^; moreover, it is invariant under the action 
of G on TQ. From the defining relations of the complete lift of a type (1, 1) tensor field above, 
one easily concludes that the horizontal distribution defined by the Vilms connection is spanned 
by the complete and vertical lifts of the horizontal vector fields of the original connection. 

Next, a remark about complete and vertical lifts. Let <j) : M ^ N he & smooth map, and suppose 
that vector fields U on M and V on N are (/^-related. Then C/^ and are T(^-related, and 
likewise C/^ and V"^ . One can easily prove this by considering the flows of the involved vector 
fields. 

Consider now (X^)^ , for any vector field X on Q/G. It is horizontal with respect to the Vilms 
connection. By the previous remark, since X^ is vr-related to X, (X^)^ is Tyr-related to X^ . 
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Thus (X^)^ is the horizontal hft with respect to the Vihns connection of the vector field 
on T{Q/G). With some abuse of notation we may write {X^)^ = {X^)^ (warning: V and H 
have different meanings either side of the equality sign). Likewise for X^. So we can rewrite 
the horizontal Euler-Lagrange equation above as follows: 

r((y^)«(L))-(y°)«(L) = o 

for all Y on Q/G. 

The Vilms connection is G-invariant, so each of (Y^)^ and (Y'^)^ is an invariant vector field on 
TQ, and so each passes to the quotient to define vector fields (Y^)^ and (Y'^)^ on TQ/G. The 
reduced horizontal Euler-Lagrange equation is 

t{{Y''f{l))-{Y''f{l) = 

on TQ/G, for all Y on Q/G. We call it the horizontal Lagrange- Poincare equation. 
In the next step, we can write this equation as 

f((y^,(i"0) - {Y'',d^i) = o. 

Here d^l is a 1-form along the projection TQ/G — ?> T{Q/G) (or a 1-form on TQ/G which 
is semi-basic with respect to that projection), such that for any vector field W on T{Q/G), 
{W,d^l) = W^{1). Since Y^ and Y"" are actually the lifts from Q/G to T{Q/G) this looks very 
much like an Euler-Lagrange equation on T[Q/G) (in fact, it would be one if d^^l was replaced 
by the exterior derivative d on T[Q/G)). 

We have divided the reduced equations in two sets, in accordance with the decomposition of 
TQ/G given by the Atiyah sequence. We conclude therefore: 

Proposition 6. The Lagrange- Poincare equations are given by 

m,p)) = {{c,p)) 

f((y^,d«0) - (y°,d*'0 = 0, 

where C is any G-equivariant Q-valued function on TQ and Y is any vector field on Q/G. 

3.2 Lagrange-Poincare-type reduction of nonholonomic systems 

Again, there is a sense in which the reduction of the Lagrange-d'Alembert equations is immedi- 
ate. These equations say that on C 

r(Z^(L)) - Z°(L) = 

for all Z T>. We assume that "D is G-invariant. It is again enough to take Z to be invariant. 
The function Z'^{L) is invariant, and so defines a function on TQ/G, which is just Z'^{1) (where 
/ is the reduced function of L). Likewise for Z^{L); however, since Z & T>, Z^ is tangent to C 
and so we can replace I by 1^ in the first term. Then the reduced equation is simply 

f (Z^(/,)) - Z^(/) = 

on C/G, for all Z € D; of course the second term must be understood as the restriction of that 
function on TQ/G to the submanifold C/G. 
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The reduced equations are now called the Lagrange-d'Alembert-Poincare equations. Like the 
Lagrange-Poincare equations they can be split into two sets, corresponding to the vertical and 
horizontal parts of T>. Much as before, these two sets are dictated by the version of the Atiyah 
sequence we found in Proposition [5j 

The Lagrange-d'Alembert-Poincare equation for momentum. First, we consider vertical 
vector fields in T>, that is, vector fields Z in S. Every such vector field defines a g-valued 
function ^ on Q, where now ^(g) € g'''. Since L is G-invariant we may define the momentum 
map p : TQ — ?> g*, as usual, and it is G-equivariant with respect to the coadjoint action on 
g*. However, we now have no reason to suppose that in general p, or any component of it, is 
conserved. On the other hand, the argument that leads to the formulae Z^{L) = {C,p) and 
Z'^{L) = {C,p) still holds good, and we conclude that the Lagrange-d'Alembert equation for 
Z ^ S can be written 

on C. We conclude further that the following weakened version of conservation of momentum 
for a constrained system holds: {(,r{p)) = for all g-valued functions C such that ({q) S g*^; 
that is to say, for all {q,u) G C, r(^,^)(p) € (s'^)^- Of course, if it should happen that for some 
C) (C)P) = on C then {C,p) will be conserved. In particular, this will occur if S contains a 
fundamental vector field of the action, that is, if there is some A € g such that A € g"^ for all 
q (z Q: then {A,p) (the ^-component of momentum) will be conserved. 

To obtain a G-invariant vector field Z G S we must take C to be G-equivariant under the adjoint 
action on g. Then we have the reduced equation 

m,p)) = ((c,p)) 

on C/G, where C{q) G g''. 

The horizontal Lagrange-d'Alembert-Poincare equation. To obtain the reduced equa- 
tion corresponding to the horizontal part of V we need a splitting of the modified Atiyah sequence 
of Proposition [5l One may derive such a splitting from a principal connection on ir : Q ^ Q/G 
with the property that the horizontal lift of V is contained in T>. Such a connection can be 
constructed locally, by defining its horizontal subspaces as follows. Take a local section of vr. 
For every q in the image of the section choose some complement to Sq in Dq and extend it to a 
complement of Vq in TqQ, smoothly over the section. Finally, extend the result along the fibres 
by the action of G. 

The reduced equation is 

f((y^,d^/)) - {Y'',(fi) = 

as before, but now with y G P. 

The conclusion of this section is therefore: 

Proposition 7. The Lagrange-d'Alembert-Poincare equations are given by 

f((C,P» = ((C,P)) 
f((y^,d«/)) - (y°,d»o = 

on C/G, where C{q) S 0'^ CLnd Y gV. 
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4 Lagrange-Poincare-type reduction: formulae 



The versions of the reduced equations given in the previous section are elegant and instructive. In 
the hterature one may find other geometric approaches to obtain the reduced equations (see e.g. 
[5tll3] for the Lagrange- Poincare equations and [ll[3l[6l[20] for the Lagrange-d'Alembert-Poincare 
equations). In the interest of comparison, we shall now formulate local versions of our reduced 
equations. Our method is entirely based on the use of suitably adapted anholonomic frames on 
Q. Unsurprisingly, the version of the Lagrange-d'Alembert-Poincare equations we finally obtain 
will combine elements of both the Lagrange-d'Alembert and the Lagrange-Poincare equations, 
so we first deal with each of those cases separately. 

But even before doing so, it is worth recalling the basic formulae relating to anholonomic frames. 
Let {Xi} be an anholonomic frame on Q and v'^ the quasi- velocities corresponding to that frame. 
(The quasi-velocities are not to be confused with the canonical fibre coordinates associated 
with the X*; the coordinates (x*,?;*) are to that extent unnatural.) A second-order field V can 
then be written in the form 

T = v'Xf + rXj. 

We write [Xi,Xj] = R^-X^, where the functions R!!ij are collectively called the object of an- 
holonomity. Let f * be the quasi-velocities corresponding to the frame: then 

Xf{v^) = -E^,,v\ Xr{y^) = 4- 

In terms of coordinates x* on Q we may write 

— yj ^ _ „.fc ^ — 

where Xi = Xjd/dx^ . The first term in the expression for X? is formally the same as Xi itself, 
but is of course a local vector field on TQ rather than a vector field on Q; we shall continue to 
denote it by Xi, though this is strictly speaking an abuse of notation. With this understood, 
the Euler-Lagrange equations may be written in Hamel form, 

r(^)-X,{L) + Rl,v'^=0 
\dv' J ^ ' dv^ 

(as in e.g. [3]), which is the form we had announced in our Introduction. 



4.1 The Lagrange-Poincare equations 

We take the frame {Xi} on Q to be invariant under the action of G, and to be of the form 
where the X^- are vertical and such that their values at any point q form a basis for 
the vertical vectors at q. Then the Xj are invariant and transverse to the fibres of Q ^ Q/G, 
and may be considered as the horizontal lifts of their projections Yj to Q/G, with respect to 
some principal connection u; note that {Yj} is a frame for Q/G, in general anholonomic. 

Let {Er} be a basis of g and let Ej. be the fundamental vector fields of the action corresponding 
to this basis. We have [Er,Es] = —G^^Et where the coefficients C*^ are the structure constants 
of Q with respect to the given basis. A vector field on Q is invariant if and only if all [X, E,.] = 0. 

The vector fields Xj. could be obtained by taking a local section of Q ^ Q/G, choosing a basis 
of vertical vectors at each point of the section varying smoothly over it, and using the G-action 
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to define the vector fields along the fibres. If one chooses the initial values of the X^- to be Er 
we obtain the vector fields Er that we have used in previous publications [9l[T9]. However, in 
view of what we need for the next section on constrained systems, it will be convenient to work 
in greater generality already in this part of the section. 

We may write Xj. = X^Eg where the coefficient matrix is nonsingular. We have 

[Er,Xs\ = [Er,XlEt] = {Er{Xl) - Cl^X^)Et, 

and so the necessary and sufficient condition for the X^ to be invariant is that the coefficients 
X* satisfy 

Er{Xs) = C^uX"^. 

It follows immediately that Xr{Xl) = C^^X^X^, whence 

where the overbar indicates the matrix inverse. Since Xr, Xg and Xt are all invariant, so must 
the coefficient be. We set C*^ = X^C^^X^X^ . Then each C^g may be treated as a function 
on Q/G, and the collection of such functions may be regarded as the structure constants of g, 
though expressed in terms of the Xr- (If Xj. = Er then (7*^ = C*^.) 

The equation EriXl) = C^^Xg expresses the fact that the g-valued function X^Et = E,g on Q 
is G-equivariant with respect to the adjoint action, and therefore corresponds to a section of 
the adjoint bundle g Q/G] thus the together form a local basis of sections of the adjoint 
bundle. Now 

Xiiir) = Xj{X:.)Eg = X!Xj{X'r)^g; 

for convenience we shall write for XfXj{X^.). We know that the functions T|,^ are G- 
invariant, and may therefore be regarded as functions on Q/G. Now Xj is the horizontal lift of 
the vector field 1/ on Q/G to Q, so we have 

TUs = Yfi^r). 

This means that the are the connection coefficients, with respect to the local basis {Yj} of 
vector fields on Q/G and the local basis {^r} of sections of g — > Q/G, of the connection induced 
by CO on the adjoint bundle (see e.g. [9] for more details). 

Since the elements of the frame {Xi} are invariant, so are their brackets, and so is each i?^-; 
it may therefore be regarded as a function on Q/G. The following facts about the R^j are 
important. 

• The constitute the object of anholonomity of the frame {Yj}. 

• The vertical component of which is R^jXr, is closely related to the curvature of 
the connection uj: in fact the curvature, as a g-valued function, is —R^jX^Er. We write 
—Kjj for R'jj as a reminder of this fact. 

• Since is always vertical, Rl^ = 0. 

• R% = '^h- 
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We let {v''' ,v^) be the the quasi- velocities corresponding to the frame {Xr,Xj}; thus for v € TqQ, 
V = v^Xr{q) + Xi{q). But then tt^,v G T^(^g-^{Q/G) is given by 

TT^v = v^-n^.Xi{q) = v^Yi{Tr{q)). 

From the invariance of the frame {Xi} we conclude that the i)* are invariant, and therefore 
constitute fibre coordinates on TQ/G Q/G. Since vr* : TqQ T^(^q-^{Q/G) is surjective, we 
can identify with Tr*w^ , where the are the quasi- velocities of the frame {Yj}. 

We now consider the Euler-Lagrange equation for momentum, T{X^{L)) — X^{L) = 0. The 
components of momentum are of course given by pr = E^{L). But we are working in terms of 
the invariant basis {^r}- Let us set Pr = X^{L) = X^ps = {£,r,p)', then Pr is invariant. The 
Euler-Lagrange equation becomes 

riPr) = X'rPs = XfXlPs, 

which is the component form of the momentum equation given in Section [31 and is usually 
referred to as the momentum equation in a moving basis [1, 2]. Now in terms of quasi- velocities 

— = v'X, = v^Xi + v''Xr, 
dt 

whence 

X^X', = Xt{v'Xi{Xl) + v'^X^iXD) = Tjy - C'rtv\ 

and we have 

r(p,,) = {Tiy - c:.,v')Ps. 

Taking account of the known facts about the R^j, together with the invariance of the frame, we 
have 



x^ 



_d_ 



d 



X 



V 



By substituting these expression in the reduced equations t{X^{l)) — X^{1) = (with i = r and 
i = I successively) we get: 

Proposition 8. The Lagrange- Poincare equations are given by 
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The first of these is of course just the reduced form of the momentum equation given earher. 

If we take the Yj to be coordinate fields the horizontal equation takes on a somewhat more 
familiar appearance: 

now the are effectively the standard fibre coordinates on T{Q/G) (quasi no longer). 

Finally, we reconcile this version of the horizontal Lagrange-Poincare equation (for any frame 
{17} on Q/G) with the more abstract one given earlier, 

r{{Y'',d^)- {Y'',d^l) = 0, 

by computing cPl in terms of the frame {Yj}. Let {1?^} be the basis of 1-forms on Q/G dual to 
the Yj, and the quasi- velocities. Then {i?^, dv^} is a basis of 1-forms on T{Q/G) (we haven't 
distinguished notationally between 1-forms on Q/G and their pullbacks to T[Q/G)). Note that 

{Yr,^'') = 5i, {Yf,dv') = -R^j,v^, (y/,^-^) = 0, {Y^,dv') = 5i. 

Now Xj = Yf, from which it follows that 



{Yf f = Xf = Yj- Rfy^ + {K}y - TL.^) ^ 



Recall that d^l is a 1-form along the projection {TQ)/G T{Q/G), which means that it may 
be expressed as a linear combination of the forms {t?^ , dv^ } with coefficients which are functions 
on {TQ)/G. Using the expressions above for (Yf)^ and (Y^)^ we obtain 

dH = (Yjil) + {K}y - T},v^) + ^dv\ 

which leads to the expressions for the reduced equations given above. 



4.2 The Lagrange-d'Alembert-Poincare equations 

We now choose the anholonomic frame {Xi} in the form {Xa,Xa} where {^a} is a local basis 
of the distribution D. We write the quasi- velocities as {v°',v"'). The constraint submanifold C is 
then simply given by v'^ = 0. The Lagrange-d'Alembert equations are T{X^{L)) — X^{L) = 
on C. Recall that T represents here a vector field on C of second-order type, which means that 
it is of the form 

r = v^x^ + rx^. 

The Lagrange-d'Alembert equations become 

in Hamel form. It is sometimes considered preferable to separate out those terms which involve 
differentiation along C from those which involve differentiation transverse to it; in the former we 
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can replace L by Lc, the constrained Lagrangian, in other words the restriction of L to C. We 
obtain 

To obtain reduced equations for an invariant constrained system we need an adapted frame 
{Xj} = {X(y,Xa} which is invariant, as before. The basis {X^} of T), in turn, is of the form 
{Xp,X^} where {Xp} is a basis for S. The set {Xa} takes the form {X^^Xk} where the Xc 
are vertical. The collection {Xp^X^} is a basis {X^} of the vertical vector fields; in general we 
can no longer take E^. for X^.. The collection {X^^X^} = {Xj} is transverse to the fibres of 
Q ~^ Q/G and is invariant, so can be taken to be the horizontal lifts of their projections Yj to 
Q/G with respect to some suitable principal connection u. The vector fields form a basis for 

v. 

The corresponding quasi-velocities are (t;",?;") or [v^ ,v'^ ,v'^,v^); the constraint submanifold 
C is given by v"" = 0, and (v'^,v^) are quasi-velocities on Q/G corresponding to the frame 
{Yj} = {Yi^,Yk}, with = defining the constraint submanifold C. 

The Lagrange-d'Alembert equations reduce, taking a = p and a = k in turn, to the following 
pair of equations on C/G: 

Now L and C are both invariant under G, and so the constrained Lagrangian is invariant 
under G, and defines a function on C/G, which coincides with the restriction of / (a function 
on [TQ)/G). The function is called the constrained reduced Lagrangian (but might just as 
well be called the reduced constrained Lagrangian). We can use this to rewrite the reduced 
equations. 

Proposition 9. The Lagrange- d'Alembert-Poincare equations are given by 




The first of these is the reduced momentum equation. The restriction of the momentum p to 
0-^, that is, the map : TQ (o^)* given by 

{C,p'^{q,u))=i^iL){q,u) for ^ G fl", 

is sometimes called the nonholonomic momentum map. Its components are the G- invariant 
functions Pp = {^p,p) = Xp{L), where ^p = XpEr defines a section of g-'^ — )• Q/G. The 
functions Pp satisfy T{Pp) = X^XpPr on C, and this reduces to the Lagrange-d'Alembert- 
Poincare equation for momentum given above. 
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We next consider some special cases of Lagrange-Poincare-type reduction of nonholonomic sys- 
tems. 



In case TqQ = Dq + Vq (i.e. when the so-called 'dimension assumption' is satisfied), the space 
C is the whole of T{Q/G). Furthermore, we can replace the Yj with coordinate fields d/dx^ on 
Q/G, and the horizontal reduced equation becomes 



d}_ 



C/G 



In case Sq = {0} the constraints are said to be purely kinematic in [7]. In this case there is no 
momentum equation. 

Chaplygin systems (see e.g. [UlIS]) are systems which have both of the above properties. There 
is no momentum equation, and P is now the horizontal distribution % ola, principal connection. 
We can therefore identify C/G with T{Q/G). The reduced vector field is now of the form 

i.e. it is a (true) second-order differential equation field on Q/G^ and its coefficients can be 
determined from the equations 



dlc_ 



dlc_ 
dx^ 



TiQ/G) 



These equations are of the form of Euler-Lagrange equations subjected to an external force of 
gyroscopic type. See e.g. [10] for more details on this case, in the framework of anholonomic 
frames. 

In case P C V there is no horizontal equation, and the momentum equation is just 

dl 



dvP J 



C/G 



One important special case occurs when the configuration space Q is a Lie group (that is, Q = G), 
and the constraints are linear; the reduced equations are then called the Euler-Poincare-Suslov 
equations in e.g. [31 E]. 



5 Routh-type reduction for systems with horizontal symmetries 

We now consider the class of systems with a so-called group of horizontal symmetries [2]. For 
that case, one assumes that there exists a subgroup H C G, the so-called group of horizontal 
symmetries, such that 1 G P for ah A G f) and Sq = Vq CiVg = Vg = {A{q) \Ae\:),qeQ}. 
Because of the property g^i^'^^ = ad(5'~^)0'^ that we encountered in Section [2l we get that 
f) = ad{g~^)[), meaning that f) is necessarily an ideal (or that H is a normal subgroup). 

For systems with the above properties, one can, of course, still use the reduction procedure as 
described in the previous sections. There are, however, also other approaches to reduction. For 
example, in [7] it is shown that a version of Marsden-Weinstein reduction can be applied to 
this case. The goal of this section is to show that one can easily stay on the 'Lagrangian side' 
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and recast everything in terms of Routh reduction. We will follow the geometric approach to 
non-Abelian Routh reduction we have developed in [8]. (For a different approach see [H].) As 
before, the main observation is that all one needs to do is to choose an appropriate frame. Let's 
assume for simplicity that in this section Vq + Vq = TqQ. 

Let {X^} be the invariant vector fields we had before. If {Er} = {Ep,Ec} is a basis of g whose 
first members {Ep} span f), we can use {-'^^q} = {Xf^,Ep} as a (now not-invariant) anholonomic 
frame for V, and {Xf^, Ep, Ec} as a complete basis of vector fields on Q (with corresponding 
quasi-velocities {v'^ , v'' , v^) . Given that Ep{L) = 0, the Lagrange-d'Alembert equation in the 
direction of Ep now becomes the conservation law 

e;{l) = ^ip, 

where n = UpE^ G [)*, where the E^ are part of the basis that is dual to {Ep,Ec}. This 
represents a relation on C, not on the whole of TQ. 

The remaining Lagrange-d'Alembert equations are of the form 

r(x:(L))-x,^(L) = o. 

We will restrict these equations to a fixed level set of momentum, from now on denoted by Np, 
and rewrite them in a form that contains only vector fields that are tangent to Np. We will do 
so in two steps. It is easy to see that the vector fields 

X^ = X^ + R:yE^ and X^ 

are tangent to C. Here R^^^ stands for the component of the bracket along Er- There 

is no contribution in v'' since [Xf^,Ep] = 0. The equations then become 

r(X:(L)) - X^^L) = -Rl,v'KiL) 

on Np. 

If we further assume the matrix (gpa) = {Ep{E^{L))) to be non-singular, the momentum equa- 
tions can be solved in the form = l^, where are functions of the other variables. Moreover, 
under that assumption we can always find vector fields and , with 

w^ = x^ + ap,e; 
w: = x: + bpe;, 

which, as well as being tangent to C, are also tangent to the level set Np (the notation may again 
be a bit misleading since they will not be complete or vertical lifts). Let us denote pp = Ep{L) 
and let us introduce the Routhian of L as the function TZ = L — v^Pp. Then, on the level set 
(which is a part of C) 

W^iTZ) = W^iL) - W^iv^Pp = X^{L) + APpp - X^iv^Pp - A^pp 
= X^{L) - X^iv^Pp = X^iL) + Kp^^f^p - Kxv'S?.f,p 
= X^,{L) + R'^^ytip - R'^^yfip = X^{L) 

W:in) = W^iL) - W^iv^Pp = Xl{L) + BP^pp - BPpp 
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Here, i?^^ stands for tlie component of [Xk,Xa] along Ep. We have also used that [X^jEp] = 0. 
The vector fields E^ are tangent to C. We can fix functions such that the vector fields 

e^ = e^^+c;e: 

are tangent to N^. 

Since F is tangent to pp = fip its restriction to this level set is of the form 

r = v^W^ + lPE^ + {V' o l)W^. 

The coefficient o /, can be determined from the remaining Lagrange-d'Alembert equations, 
which take the form 

nwnnn) - wnnn = -Rixv^E^m - <,^Vp 

on iV^. 

We can now try to understand how to reduce this restriction of F. It is easy to see that the 
action of G on C can be restricted to an action of the isotropy group Hp^ on the level set Np in 
C. Indeed, we have 

= A-E'^iE^iL)) = -A-ClpE^^iL) = -A^Clpiir 

if and only if ^ = A'^E^ € i)p. We can therefore reduce the above vector field to a vector field 
Ti on Np/Hp. This reduction method is the direct analogue of the situation for standard Routh 
reduction (in the absence of constraints). 

But there is more. Since we know that H is normal in G, the level set Ep{L) = fip has also the 
following behaviour 

= A^E^{e;{L)) = -A^C^E^iL) = -A^-C^^E^iL) = -A^C^^f,^ 

if and only A = A^E,. € 0^. Therefore, the G-action on C restricts in fact to a G^-action on 
the level set Np. We are now in the situation of a G- invariant vector field F on a manifold C, 
which we can restrict to a G^-invariant vector field on Np and which we can reduce to a vector 
field fa on Np/Gp. 

The link with the vector field Fi of the previous paragraph is the following. Instead of doing a 
direct reduction by Gp, one can perform a reduction in two stages. Indeed, it is easy to define an 
action of Gp/Hp on Np/Hp (see also [7]). The vector field Fi will be invariant under that action 
and we can therefore perform a second reduction by means of its symmetry group Gp/Hp. 

We will not write down explicit expressions for these reduced vector fields and their corre- 
sponding differential equations. Instead, we will make the situation clear by means of a simple 
example. 

Example. Consider the system with L = ^(x^ + + i^) on M'^ with constraint z = xx 
(a variation on the theme of a nonholonomic particle). This example is taken from [7], but 
we will rephrase it in our current framework. We have V = span{5/5x + xd/dz,d/dy} and, 
since the system is invariant under the R^-action given by (r, s) x (x, y, z) i— > {x,y + r, z + s), 
V = span{d/dy,d/dz}. Therefore, S = T> (1 V = span{d/dy}. This coincides with the case 
where H = M x {0}, with action (r, 0) x (x,y,z) i— )■ (x,y + r,z). Remark that V + V = TQ. 
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Quasi-velocities with respect to the given frame are Vx = x, Vy = y and Vz = z — xx. The vector 
field = X = d / dx + xd/dz is invariant under the G-action. 

The preserved momentum is here 

The remaining equation on C is T{X^{L)) - X'^{L) = 0. Since ^^(u^) = and X'^{vz) = 0, 
both X^ and X'-^ are tangent to the constraint, so we can rewrite that equation as r(X^(Lc)) — 
X'^{Lc) = 0, with Lc = i((l + x^)iP' + ij^). One easily verifies that this equation is equivalent 
with 

[1 + x )x — xx =0. 
This equation is evidently M^-invariant and the reduced vector field is 

^ . d xx^ d 

dx 1 + x^ dx 

It is instructive to see how one gets the same result when we use Routh reduction. Remark that 
if^ = M X and = M2. Since also = and ^^(y) = the vector fields = X^ 

and X'^ = X'^' are already tangent to the level set ij = ji. We can therefore simply re-write the 
remaining equation as 

^(x-(7^p)-x°(7^;i) = o, 

where TIP^ is the restriction of the Routhian IZ^ to the constraints and to the level set. It is 
given by 

7^^ = i((l + xV-/^')• 

The Routh equation above is again (1 + x^)x — xxP' = 0. We can now do a direct reduction by 
means of the largest group G^. We see that N^/G^ = TM. Due to the absence of gyroscopic-type 
terms, the G^-reduced version of the above equation will be a genuine Euler-Lagrange equation, 
with the G^-reduced Routhian as its Lagrangian. The vector fields X^ and X'^ reduce to the 
vector fields d/dx and d/dx on M and the reduction actually amounts to cancelling the cyclic 
variables y and z from the above equation. A similar reasoning holds for the reduction in two 
steps. We have = R and G^/H^ = M; the first reduction cancels y and the second cancels z. 
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